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Galileons are higher-derivative theories of a real scalar which nevertheless admit 
second order equations of motion. They have interesting applications as dark en- 
ergy models and in early universe cosmology, and have been conjectured to arise 
as descriptions of brane dynamics in string theory. In the present paper, we study 
the bosonic sector of globally = 1 supersymmetric extensions of the "standard" 
Galileon Lagrangians in detail. Supersymmetry requires that the Galileon scalar 
now becomes paired with a second real scalar field. We prove that the presence 
of this second scalar causes the equations of motion to become higher than second 
order, thus leading to the appearance of ghosts. Our result suggests that super- 
symmetric Galileons should be regarded essentially on the same footing as other 
higher-derivative terms arising in effective field theories. We also analyze the energy 
scales up to which, in an eff'ective description, the ghosts can be tamed. 
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I. INTRODUCTION 



Galileon theories of a real scalar field are special because they have two-derivative equa- 
tions of motion despite having higher-derivative Lagrangians. They are a sub-class of the 
most general scalar theories with two-derivative equations of motion, known as Horndeski's 
theories [ll (see also fl). The "standard" Galileons |3[ have the additional property that 
in the equations of motion there are precisely two derivatives acting on each field. An im- 
mediate consequence is that the standard Galileons are invariant under a so-called Galilean 
shift symmetry </>—;■ + c + b^x^ with c, 6^ being constants, whence they derive their name. 
The property of having equations of motion with no more than two derivatives acting on 



a field is crucial since it helps to evade Ostrogradsky's theorem 



J]. That is, despite the 



higher- derivative nature of the Lagrangians, these theories do not contain ghosts. 

Galileons have attracted considerable interest due to their rather remarkable properties. 
For example, they admit de-Sitter-like solutions in the absence of a cosmological constant 
Q] and they lead to a Vainshtein-type screening mechanism so that they can be in agreement 
with solar system "fifth force" constraints while contributing a fifth force on large scales j^,^]. 
Moreover, they allow for solutions that violate the null energy condition without leading 
to the appearance of ghosts 10|, This last property means that Galileons also have 



applications to early universe cosmology, allowing the construction of emergent cosmologies 
(see, for example, the model of Galilean genesis 1^) and non-singular bouncing cosmologies 
such as new ekpyrotic theory 13l-ll8| or the matter bounce model 191. Such alternative 



models to inflation even play a significant role in eternal inflation 



20 



22|. 



There exists a suggestive construction of Gali 
ing the dynamics of co-dimension one branes 



eon Lagrangians as the theories describ- 



23|. This has led people to speculate that 



Galileons might arise naturally out of string theory and, hence, enjoy a more fundamental 
status than other higher- derivative terms, in analogy to the Dirac-Born-Infeld action. Brane 
backgrounds in string theory typically preserve some amount of unbroken supersymmetry. 
Therefore, if Galileons are to arise from string theory it will be in a supersymmetric context. 
Hence, it is of importance to study the supersymmetric extensions of Galileon theories. In 
previous work 2J], it was shown that conformal Galileons can be made globally = 1 
supersymmetric-these theories arising naturally as a way of obtaining correct sign spatial 
gradients in supersymmetric ghost-condensates (see also 25|, l26|). It was found that the 
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new fields required by supersymmetry (a second real scalar, a spin ^ fermion and a com- 
plex auxiliary field) admit stable, positive-energy fluctuations around specific backgrounds, 
namely those where the second scalar field is constant. However, possible ghost instabili- 
ties associated with vacua with a spacetime dependent second scalar were not explored. We 
will do this in the present paper, restricting our discussion, for simplicity, to the standard 
Galileons of js] within the context of four- dimensional global = 1 supersymmetry. 

To begin, we present complex scalar Galileons which, when the second scalar is set to zero, 
reduce to real Galileons of the L3, L4 and L5 type. These possess manifestly two-derivative 
equations of motion and a Galilean symmetry for the two constituent real scalars fields. 
We then show, however, that such complex Galileon theories cannot be obtained in = 1 
supersymmetry. We next consider the cubic-in-the-field, four-derivative L3 Lagrangian, and 
show that there is a unique possible = 1 supersymmetric generalization. However, it is 
demonstrated that this Lagrangian leads to higher- derivative equations of motion] An im- 
mediate consequence is that, around general backgrounds, this theory admits a ghost, whose 
existence we explicitly demonstrate. In the effective field theory context, we then calculate 
the mass of the ghost and argue that for a sufficiently low cut-off scale the ghost degree of 
freedom can be safely ignored. In the final technical section before the discussion, we ex- 
tend our analysis to supersymmetrize the quartic-in-the-field, six-derivative L4 Lagrangian. 
Again, it is found to lead to higher-derivative equations of motion. This time we perform 
the stability analysis using the canonical Hamiltonian formalism, and explicitly demonstrate 
the existence of ghosts as well as the unboundedness of the Hamiltonian. Our results imme- 
diately generalize to supersymmetrized L5 as well. The implication is that supersymmetric 
Galileons cannot be considered as fundamental, but must be treated in an effective field 
theory context in the same manner as generic higher-derivative terms. It is notable that, for 
once, the inclusion of supersymmetry does not improve the stability properties of a theory 
- quite to the contrary! 

We note that our conclusions only apply to L^, L4 and L5, but not to Li,L2 since the 
latter two Lagrangians do not contain higher derivatives to start with. Furthermore, we have 
performed our analysis within the context of global rather than local supersymmetry. How- 
ever, the generic supersymmetric structure of the higher-derivative scalar field Lagrangians 



is not substantially altered in the presence of gravity (see e.g. 27|, l28|). That is, the exis- 



tence of ghosts in the L3, L4 and L5 Galileons will persist when these are coupled to A^ = 1 
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supergravity. Finally, recall that our results are derived for the "standard" Galileon theories. 

io|, 



Many variants of the original model have been constructed, such as conformal Galileons 

, , bi-Galileons 



DBI Galileons 
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32| and so 



Galileons with an internal symmetry 
on. We do not have a proof that our results extend to these variants as well, but preliminary 
calculations strongly suggest that they do. This will be discussed elsewhere. 



In this and the fol 
Lagrangian given by 



II. GALILEONS AND COMPLEX FIELDS 

owing two sections, we will focus on the simplest non-trivial Galileon 

^3 = -\{d<Pf^(t>- (II.l) 



By varying with respect to 0, one can immediately see that the equation of motion is second 
order and given by 

(□0)2 _ = 0. (II.2) 

Thus, despite the higher- derivative nature of the Lagrangian, the equation of motion is well- 
behaved and the Cauchy problem is well-posed. In four dimensions, there are two more such 
Galileon Lagrangians, 

U = -^(90)2((n0)2-0''^>,^,), (II.3) 

L5 = -^{d<pf{in(t)f - 3n#''^'^0,^, + 20'^>,^,0,/) (II.4) 

which also lead to second-order equations of motion. For example, the equation of motion 
for L4 is given by 

(□0)3 - 3n00'^^0,^, + 20''^>,^>,/ = 0. (II.5) 

A detailed discussion of the L4 Lagrangian and its supersymmetric extension will be dis- 
cussed below in the final technical section. 

In = 1 supersymmetry, scalar field theories can be constructed using chiral superfields 
$. The lowest component of such a superfield is a complex scalar A, which can be decomposed 
into two real scalars as 

A=^{<P + iO. (II.6) 
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One consequence is that supersymmetric scalar field actions can always be written as her- 
mitian combinations of A and its complex conjugate A*. At this point, it is interesting to 
note that there is an immediate generalization of Galileons to complex Galileons. These are 
obtained by replacing (p — )■ \/2A and then taking the real part. For L3 above, this amounts 
to considering the Lagrangian 

L^ = — ^^{dAfDA + h.c. , (11.7) 
v2 

where h.c. stands for "hermitian conjugate". It is then evident that the resulting equations 
of motion are still second order, since they are given by 

{DAf - A'^^A^^^ = 0, {nA*f - A*''"'A*^^ = 0. (II.8) 

In terms of the real scalars (p and ^, the Lagrangian and equations of motion are 

^3 = -l{id<P)'a<p - idO'a<p - 2d<p ■ d^n^) , (11.9) 

= (□0)2 - - (DO' + ^'''i.u, (11.10) 

= - ^•'^'^Cm-, (11.11) 

clearly exhibiting that we now have a coupled two-field Galileon system. Not only are the 
equations of motion of second order, but both fields admit independent Galileon-type shift 
symmetries (p ^ (f) + c'^^^ + b^/f^x'^ and — )■ ^ + c^^^ + b^^x'^ respectively. 

Note that other actions, involving both A and A* in a single term, do not lead to second- 
order equations of motion. To illustrate this important point, consider the action 

= -^dA-dA*DA + h.c. (11.12) 

= ~{id<pfn<P+idO'a<p), (11.13) 

leading to the equations of motion 

= (00)2 - - e^^^Mv - u/', (11.14) 

= D^D^ + (11.15) 
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Clearly, these are higher-order in time and, thus, by Ostrogradsky's theorem {4], lead to the 
appearance of ghosts. 



III. SUPERSYMMETRIC GALILEONS 



In this section, we will construct all possible supersymmetric Lagrangians involving the 
product of three fields and four space-time derivatives, in order to see if there might exist 
inequivalent supersymmetric extensions of the L3 Lagrangian (lll.ll) . We will work in = 1 
superspace (for a detailed exposition see jssi). Here, in addition to ordinary four- dimensional 
bosonic spacetime one adds four fermionic, Grassmann- valued dimensions. These have coor- 
dinates da and transforming as a two-component Weyl spinor and conjugate Weyl spinor 
respectively. One can then define the superspace derivatives 

which satisfy the supersymmetry algebra 

{D,,Dj = -2i<^a^ . (III.2) 

Any superfield can be expanded in the ant i- commuting coordinates 6', 9, with the expansion 
terminating at order 9999 because of the Grassmann nature of the fermionic coordinates. A 
chiral superfield $ is defined by the constraint 

D<l> = . (111.3) 



This has the expansion 

$ = A{x) + \/29x{x) + eeF{x) 
+i9a"'9dmA{x) - -^99dmX{x)cr"'9 + -9999UA{x) , (111.4) 

where A is a complex scalar, Xa is a spin-| fermion and F is a complex auxiliary field. In 
this paper, we will ignore the fermion. Furthermore, since we are only interested in the 
structure of kinetic energy terms, we need not introduce a superpotential - in the absence 
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of which the F field can, and will, be consistently set to zero. 

What makes superspace so useful is that the top component (that is, the 6666 com- 
ponent) of a superfield transforms under supersymmetry into a total spacetime derivative. 
Hence, one can use this top component to construct supersymmetric Lagrangians. The 
top component can be isolated by integrating the superfield Lagrangian over superspace 
with d'^6 = d'^dd'^d or, alternatively, by acting on it with D^D^. The supersymmetry algebra 
(1III.2|) then implies that the top component of a superfield will contain two additional space- 



time derivatives compared to its lowest component or compared to the superfield expression 
itself. For example, ordinary two- derivative scalar field theories are obtained by isolating 
the top component of the Kahler potential, which is an hermitian function of the chiral 
superfield $ and its hermitian conjugate $^ involving no spacetime derivatives. 

In our case, we are interested in Lagrangians involving the cubic product of a scalar field 
and four spacetime derivatives. This means that we should consider all possible superfield 
expressions involving the cubic product of a chiral superfield and two spacetime derivatives 
(and linear combinations of all such terms). The superfield Lagrangians of potential interest 
are straightforward to write down. They are given by the 6666 components of the following 
expressions (where derivatives act only on the immediately following superfield): 

df'^d^^^ + /i.e. (III.5) 
a/^$a^$t$ + (III.6) 
aM$a^$$t + (III.7) 

All other terms of potential interest can be related to these using integration by parts. 

Note that only the first of these expressions can possibly lead to the complex Galileon 
Lg given in (1II.7P of the previous section. This follows from the fact that it is the sole term 
containing only $'s or only $^'s in a single term. However, the chirality of $ immediately 
implies that this term is, in fact, zero. To see this, instead of integrating over d^^, one can 
make use of the Grassmann nature of the 6, 6 coordinates and replace d^d by a D^D^ deriva- 
tive of the corresponding superfield expression. Since D commutes with partial derivatives, 
it immediately follows that superfield expressions constructed exclusively out of $'s and 
partial derivatives must vanish, since the D derivative will necessarily act on a chiral field 
$ thus yielding zero. Note that this no-go argument relies solely on holomorphicity and. 
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thus, also extends to potential super symmetric extensions of complex Galileons with higher 
powers of fields, such as L4 and Lf. Hence, it is now clear that = 1 supersymmetry 
forbids a complex extension of the L3 Galileon with two-derivative equations of motion! 
Instead, supersymmetry yields the following actions (and their hermitian conjugates): 

y dWe9''<l>a^$t$ = J d^x{- ADADA* -DA*{dAY), (III.8) 
J dSd^^a'^^a^^^t = J d^xDA*{dAf. (111.9) 

In presenting these component actions, we have used integration by parts to simplify them 
as much as possible. The first term on the right-hand side of ( 1111.8^ is uninteresting for our 
present purposes since, even for the first real scalar field (p, it does not lead to a Galileon 
Lagrangian. Hence, we are left with a single possible supersymmetric extension of the 
L3 Galileon Lagrangian, namely the real part of ( 1IIL9I) . We note that this Lagrangian is 
equivalent to the supersymmetric Galileon Lagrangian used in 2^- Thus, we define the 
supersymmetric extension of L3 as 



^susY ^ __ J d4 6»a^$a^$$t + h.c. 

= — ^^nA*{dAf + h.c. 
v2 



1 

7 
1 

= -^{{dcpfn(i)-{dO'n(i) + 2d(i)-d^n^) . (in.io) 

Compared to the complex Galileon (111.9^ . only the sign of the last term has changed! Nev- 
ertheless, this has profound consequences, since the resulting equations of motion are now 
of third order in derivatives. They read 

= (□0)2 - + (DO' + e^^^Mv + SCmC/^ (III.ll) 

= C''"'cP,,. + i,(l>/''. (III.12) 



As one can clearly see, it is the presence of the second scalar ^ that induces the dangerous 
higher-derivative terms. We will show explicitly in the next section that the presence of 
these higher derivatives leads to the appearance of a ghost. 



9 



IV. HIDING FROM THE GHOST 

We would now like to explicitly demonstrate the ghost degree of freedom in Lf^"^^. The 
presence of a ghost is already implied by Ostrogradsky's theorem ^ and we will, in fact, 
analyze Lf ^"^^ from this point of view in the following section. Nevertheless, we prefer to 
also analyze the Lagrangian Lf ^'^^ directly, both because it is instructive to see the ghost 
appearing at the level of the Lagrangian and because such an analysis elucidates in what 
regime the ghost can be harmless. For this purpose, it suffices to look at the time-derivative 
terms in the Lagrangian, since it is these that are associated with ghosts. Adding a canonical 
kinetic term Lf^"^^ = J d'^^$$^ = —d'^Adfj^A*, as well as an overall constant C3 in front of 
the L^^^^ Lagrangian, the Lagrangian of interest becomes 

where we have integrated by parts in order to place all double derivatives on (p rather than 
^. Note that this is a completely arbitrary choice and does not reduce the generality of our 
analysis. We consider a time-dependent background and would like to study perturbations 
around it. Thus, we define 

= ^{t) + 6(P{x^) , e = m + sa^n ■ (iv.2) 

Even though the perturbations depend on both time and space, we will only be interested in 
the time dependence here. To quadratic order in fiuctuations, the Lagrangian then becomes 

L^'ln.^ = liS'<P)' + ^(1 + 2c3|)(50' + 2c^id^6'<p. (IV.3) 
By defining a new fiuctuation variable 

5b = 6^ + ^^^^ J (j) , (IV.4) 

1 + 2C30 
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the quadratic Lagrangian can then be diagonahzed to become 

^ ^ (1 + 2C30)^ 

Note that {Sty and (50)^ enter with opposite signs and, hence, one of these two terms is 
ghost-hke^. Assuming that the factor (1 + 2c3(j)) is positive, the ghost then resides in 6(f). 
As the Lagrangian shows, the significance of the ghost is essentially controlled by the size 
of C3^. This can be confirmed by looking at the dispersion relation of 6(1). If one denotes the 
four-momentum of 6(f) by p^, then the associated dispersion relation is given by 

pI{1--^^pI)=0, (IV.6) 

(1 + 2C3(^) 

where we have assumed that C, and (f> are slowly varying. The mass m is defined via p"^ = 
— Pq = —m^ and, hence, the dispersion relation implies that 6(f) consists of two modes. The 
first is a massless mode which arises from the ordinary correct-sign kinetic term. The second 
is the ghost, which has a mass 

= (i±^. (,V.7) 

44e 

Note that, as there is an overall wrong sign for the ghost in the Lagrangian, this mass is 
formally tachyonic. However, it is important to realize that this mass does not arise from a 
potential, but rather from the kinetic term {6(f))^. The implication is that this mass does not 
indicate a time scale over which the (perturbative) vacuum becomes unstable, but rather 
an energy scale associated with the ghost. In other words, as long as we are considering 
fluctuations with energy below rrig, the ghost does not get excited. From an effective field 
theory point of view, we are protected from the catastrophic instabilities associated with the 
ghost if we take the cut-off A of the effective field theory to lie below rrig. At the same time, 
we must ensure that the background itself, that is, ^, remains within the range of validity 
of the effective theory. Hence, an additional requirement is that |.^| < A^, and similar 
inequalities must also hold for higher time derivatives of ^. Together with the requirement 



^ This ghost was not seen in [2J| because in that paper the perturbation analysis was performed solely 



around ^ — constant backgrounds. 
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A < nig, this implies that we must impose (assuming {c^cpl ^ 1) 

iei<^, ili<^> ••• (IV.8) 

in order to safely suppress the ghost. Thus, as expected, for general backgrounds one must 
take the prefactor of the Galileon term to be small for consistency. 

V. L4 AND THE HAMILTONIAN ANALYSIS 

In this section, we discuss in detail the four-field Galileon Lagrangian I/4 and its super- 
symmetric extension. This Lagrangian was presented in Eq. (111.31) for a single real scalar (p. 
First, we want to point out that the holomorphicity argument presented in section UTTl im- 
plies that it is impossible to construct a corresponding supersymmetric version with purely 
two-derivative equations of motion for both real scalars 0, ^. However, it is not hard to find 
a supersymmetric extension of L4 after rewriting it using integration by parts: 

= -]d>'id<)))%id<py + ldyd^idci))''n(i) . (v.i) 

4 z 



Making use of the "building blocks" 24 



= -midAf, 1)2$ = -mOA, (V.2) 
it is straightforward to write a supersymmetric extension of ( IV. II) given by 

j^susY ^ f d4^C-^aA'(D<|.D$)a^(D$tD<|.t) + —d''^df,{D^D(^)D^(S>A + /i.c.(V.3) 

= -d^'idAydf.idA*)^ + d^Adf,{dAynA* + d^'A*df,{dA*ynA . (V.4) 

As in the previous Lf^^^ analysis, we will add this term-now with a coefficient C4-to the 
standard kinetic term Lf'^'^^ = J d^6'$$"'' = —d^Adfj^A*. Since we are primarily interested 
in the issue of ghosts, we need only consider the time-dependent part of the resulting La- 
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grangian. This is given by 

Ll^r = ^f'^'' + C4Lf (V.5) 
= -d^Ad^A* + C4( - d^'idAfdf.idA*)^ + d^'Ad^idAfDA* + d^'A*d^{dA*fDA) 

= + + 2c4C'e' + 2c4^20l (V.6) 

Again, the higher-derivative nature of the Lagrangian is manifest. The fourth-order equa- 
tions of motion that follow from this Lagrangian are 

= -■<p + Ac,^^{e'<P), (V.7) 

= -I + 4c4^(e'o - 4c4^(ee' + e^')- (v.s) 



It is instructive to carry out a Hamiltonian analysis of this theory following [3J]. For our 
canonical coordinates, we will choose ^, ^ as well as a = 4> and b = ^. The corresponding 
momenta are 

50 df dt'^ 

dL2A d dL2A I . ^2 '^^ ^ A ^'^1 . A ^72 /t r 1 ^\ 

TT. = ^ - ^ = ,^ - 4C4rTT<^ - 4C4.^r + 4C4<^0^ (V.IO) 

vr. = ^ = 4c4^^^, (V.ll) 
vr. ^ ^ = 4c4e^|. (V.12) 

The Hamiltonian is given by 

H = ^n^ + + ana + f'TTj, - Lf^f^, (V.13) 

which can be re-expressed in terms of the canonical coordinates and momenta as 

H = an, + K + ^i^r + - la^ - h\ (V.14) 

8c4 8c4 2 2 

Note that this expression is regular at 6 = since TTa and iTh both contain factors of 6^. 
To check the consistency of this analysis, one should verify that the Hamilton evolution 
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equations (j) = ^ = . . . and tt^ = — tt^ = — . . . lead to sensible results. In 
fact, the evolution equations for the coordinates are easily seen to be satisfied. Those for 
the TT^ and tt^ momenta result in 

TT,^ = 0, TT^ = . (V.15) 

These two equations are equivalent to the Euler-Lagrange equations of motion flV.7l) and 
(IV.Sp . The two remaining equations are 



7ra = a-7r0, ^^b = h - + ^-^{nl + ttD, (V.16) 



which are equivalent to the definitions of the momenta tt^ and tt^ given in f IV.QP and (IV.lOp . 



Thus, we may trust our derivation of the Hamiltonian f lV.141) . Crucially, the Hamiltonian 
depends linearly on both -k^ and vr^ and, therefore, can be made arbitrarily positive or 
negative by choosing appropriately large an^,b7r^ terms. Thus, the energy is unbounded 
from both above and below. This is a clear indication that this theory, taken literally, does 
not admit any vacuum at all and is, thus, unphysical. This explicitly demonstrates the 
presence of a ghost degree of freedom (here, in fact, there are two ghosts), and leads to 
conclusions similar to those of the Lagrangian analysis performed in Section IIV[ 

A final comment: the higher-derivative terms that we have discussed all arise because of 
the presence of the second real scalar field C,- If we momentarily fix ^ = b = 0, then the 
Hamiltonian reduces to the simple form 

%o = (v.ir) 

which is manifestly positive. Thus, by shutting one's eyes to the presence of the second field, 
one may-mistakenly-think that these theories admit a stable vacuum. Even though this 
restriction leads to fallacious conclusions when treating the above theory on a fundamental 
level, it nevertheless supports the conclusion that from an effective field theory point of view 
perturbations of sufficiently low energy around ^ = backgrounds can be admissible. 
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VI. DISCUSSION 

The fact that = 1 supersymmetric standard Gahleons containing the product of three, 
four and five chiral fields necessarily admit higher-derivative equations of motion implies that 
these theories contain ghosts. This means that when supersymmetry is included, Galileons 
lose their special status among higher-derivative scalar theories and should be treated in 
much the same way as other higher- derivative terms. That is to say, they should be re- 
garded as correction terms in a perturbative, effective field theory framework. Regarding 
supersymmetric Galileons as fundamental is untenable since their Hamiltonian is unbounded 
from above and below. It follows that these theories do not admit a stable vacuum. Although 
we do not have a generic proof, based on individual examples we strongly suspect that our 
results carry over to the = 1 supersymmetrization of generalizations of the Galileons, 
such as the conformal and DBI Galileons. By extension, they are also likely to apply to the 
non-trivial parts of Horndeski's most general scalar-tensor theory jl|. We stress that our 
results are in the context of minimal A^ = 1 supersymmetry. It would be interesting to carry 
out a similar analysis for extended supersymmetries. 

As discussed in the introduction, the brane construction of Galileon Lagrangians sug- 
gested that they could arise as the sole constitutents of membrane worldvolume theories 
in string theory-that is, in a well-defined ultraviolet framework. However, when explicit 
calculations of higher-order corrections to brane dynamics were carried out-in the non- 



M-theory sjlll -it was found that, in addition to the Galileon terms, other higher-derivative 
terms occur. These new terms are not naturally suppressed relative to the Galileons and 
lead to higher-order equations of motion. This paper shows that, with hindsight, this result 
is unsurprising-since in a full supersymmetric context the Galileon terms themselves already 
admit higher-derivative equations. 

In a fundamental framework, the presence of ghosts is catastrophic. However, in a per- 
turbative framework ghosts can be avoided, provided one looks only at fluctuations below 
a certain energy scale. An open question raised by the present work is then whether or not 
the attractive properties of certain solutions to the Galileon theories-such as the Vainshtein 
mechanism and consistent violations of the null energy condition-can be maintained in a 
supersymmetric perturbative context. We leave this question for future work. 



supersymmetric case of AdS space 35[ and in the A^ = 1 supersymmetric context of heterotic 



15 



Acknowledgments 

We would like to thank Justin Khoury and Mark Trodden for sharing their wide-ranging 
knowledge of higher-derivative theories with us. M.K. and J.L.L. would like to thank the 
University of Pennsylvania for hospitality while this work was undertaken. M.K. and J.L.L. 
gratefully acknowledge the support of the European Research Council via the Starting Grant 
numbered 256994. B.A.O. is supported in part by the DOE under contract No. DE-AC02- 
76-ER-03071 and the NSF under grant No. 1001296. 



[1] G. W. Horndeski, Int.J.Theor.Phys. 10, 363 (1974). 

[2] C. Deffayet, X. Gao, D. Steer, and G. Zahariade, Phys.Rev. D84, 064039 (2011), [1103.3260]. 

[3] A. Nicolis, R. Rattazzi, and E. Trincherini, Phys.Rev. D79, 064036 (2009), [0811.2197]. 

[4] M. Ostrogradsky, Mem. Ac. St. Petersbourg IV 4, 385 (1850). 

[5] F. P. Silva and K. Koyama, Phys.Rev. D80, 121301 (2009), [0909.4538]. 

[6] A. De Felice and S. Tsujikawa, Phys.Rev.Lett. 105, 111301 (2010), [1007.2700]. 

[7] N. Chow and J. Khoury Phys.Rev. D80, 024037 (2009), [0905.1325]. 

[8] A. Vainshtein, Phys.Lett. B39, 393 (1972). 

[9] C. Deffayet, G. Dvali, G. Gabadadze, and A. L Vainshtein, Phys.Rev. D65, 044026 (2002), 
[hep-th/0106001]. 

[10] A. Nicolis, R. Rattazzi, and E. Trincherini, JHEP 1005, 095 (2010), [0912.4258]. 

[11] K. Hinterbichler, A. Joyce, J. Khoury and G. E. Miller (2012), [1212.3607]. 

[12] P. Creminelli, A. Nicolis, and E. Trincherini, JCAP 1011, 021 (2010), [1007.0027]. 

[13] E. L Buchbinder, J. Khoury and B. A. Ovrut, Phys.Rev. D76, 123503 (2007), 
[hep-th/0702154]. 

[14] E. L Buchbinder, J. Khoury and B. A. Ovrut, JHEP 0711, 076 (2007), [0706.3903]. 

[15] E. L Buchbinder, J. Khoury and B. A. Ovrut, Phys.Rev.Lett. 100, 171302 (2008), 
[0710.5172]. 

[16] J.-L. Lehners and S. Renaux-Petel, Phys.Rev. D80, 063503 (2009), [0906.0530]. 

[17] J.-L. Lehners, Adv.Astron. 2010, 903907 (2010), [1001.3125]. 

[18] J.-L. Lehners, Class. Quant. Grav. 28, 204004 (2011), [1106.0172]. 



16 



[19] Y.-F. Cai, D. A. Easson, and R. Brandenberger, JCAP 1208, 020 (2012), [1206.2382]. 

[20] J.-L. Lehners, Phys.Rev. D86, 043518 (2012), [1206.1081]. 

[21] M. C. Johnson and J.-L. Lehners, Phys.Rev. D85, 103509 (2012), [1112.3360]. 

[22] Z.-G. Liu and Y.-S. Piao (2013), [1301.6833]. 

[23] C. de Rham and A. J. Tolley, JCAP 1005, 015 (2010), [1003.5917]. 

[24] J. Khoury, J.-L. Lehners, and B. A. Ovrut, Phys.Rev. D84, 043521 (2011), [1103.0003]. 

[25] J. Khoury, J.-L. Lehners, and B. Ovrut, Phys.Rev. D83, 125031 (2011), [1012.3748]. 

[26] M. Koehn, J.-L. Lehners, and B. Ovrut (2012), [1212.2185]. 

[27] M. Koehn, J.-L. Lehners, and B. A. Ovrut, Phys.Rev. D86, 085019 (2012), [1207.3798]. 

[28] M. Koehn, J.-L. Lehners, and B. A. Ovrut, Phys.Rev. D86, 123510 (2012), [1208.0752]. 

[29] K. Hinterbichler, M. Trodden, and D. Wesley, Phys.Rev. D82, 124018 (2010), [1008.1305]. 

[30] G. Goon, K. Hinterbichler, and M. Trodden, JCAP 1107, 017 (2011), [1103.5745]. 

[31] C. Deffayet, S. Deser, and G. Esposito-Farese, Phys.Rev. D82, 061501 (2010), [1007.5278]. 

[32] A. Padilla, P. M. Baffin, and S.-Y. Zhou, JHEP 1012, 031 (2010), [1007.5424]. 

[33] J. Wess and J. Bagger (1992). 

[34] R. P. Woodard, Lect.Notes Phys. 720, 403 (2007), [astro-ph/0601672]. 

[35] J. Khoury, B. A. Ovrut, and J. Stokes, JHEP 1208, 015 (2012), [1203.4562]. 

[36] B. A. Ovrut and J. Stokes, JHEP 1209, 065 (2012), [1205.4236]. 

[37] A. Lukas, B. A. Ovrut, K. Stelle, and D. Waldram, Phys.Rev. D59, 086001 (1999), 
[hep-th/9803235]. 

[38] A. Lukas, B. A. Ovrut, K. Stelle, and D. Waldram, Nucl.Phys. B552, 246 (1999), 
[hep-th/9806051]. 

[39] A. Lukas, B. A. Ovrut, and D. Waldram, Phys.Rev. D59, 106005 (1999), [hep-th/9808101]. 

[40] A. Lukas, B. A. Ovrut, and D. Waldram, JHEP 9904, 009 (1999), [hep-th/9901017]. 

[41] R. Donagi, B. A. Ovrut, and D. Waldram, JHEP 9911, 030 (1999), [hep-th/9904054]. 



